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Abstract 

A discrete 4>^ system is proposed which preserves the topological lower bound on the kink en- 
ergy. Existence of static kink solutions saturating this lower bound and occupying any position 
relative to the lattice is proved. Consequently, kinks of the model experience no Peierls-Nabarro 
barrier, and can move freely through the lattice without being pinned. Numerical simulations 
reveal that kink dynamics in this system is significantly less dissipative than that of the con- 
ventional discrete 0* system, so that even on extremely coarse lattices the kink behaves much 
like its continuum counterpart. It is argued, therefore, that this is a natural discretization for 
the purpose of numerically studying soliton dynamics in the continuum 0^ model. 

1 Introduction 

Solitons have diverse applications in many branches of physics. Those with applications in high 
energy physics are topological solitons which arise almost exclusively in non-integrable Lagrangian 
field theories (integrability appears incompatible with Lorentz invariance in spacetimes of realistic 
dimension) . In the absence of integrability, there is no hope of solving the multisoliton initial value 
problem exactly. One approach to the study of soliton dynamics is to perform numerical simulations 
of the field equations. In so doing one is forced to discretize space in some way, and this inevitably 
introduces fictitious discretization effects into the dynamics, which one should seek to minimize. 

The standard discretization of a field theory is obtained by replacing spatial partial derivatives 
by simple difference operators, leaving it otherwise unchanged, so that one is really studying a 
large network of oscillators, each moving in an identical substrate potential, with nearest neighbour 
couplings. Such systems are intrinsically interesting (they have applications in condensed matter 
and biophysics) and have been extensively studied in recent years |Q, ^, ^. One discretization effect 
is that static solitons can no longer be centred on any point in space: they must lie in the centre 
of a lattice cell, or exactly on a lattice site, and the energies of these two types of static solution 
are different. So there is an energy barrier, called the Peierls-Nabarro (PN) potential resisting the 
propagation of a soliton from cell to cell. As a soliton moves through the lattice, its periodic motion 
in and out of the PN potential excites small amplitude travelling waves (radiation) which travel in 
its wake, draining its kinetic energy. The soliton decelerates under this dissipation until it no longer 
has sufficient energy to surmount the PN barrier. It is then pinned to a lattice cell. 

If one is to use the standard discretization to simulate continuum soliton dynamics, these extra- 
neous discretization effects must be rendered insignificant by using a very fine spatial lattice, with 
spacing h of order (soliton width)/20. This works well, but is computationally expensive, especially 
in high dimensions. In this paper we consider an alternative approach, originated by Richard Ward, 
in the case of 0^ theory in 1 -I- 1 dimensions. Exploiting the non-uniqueness of the discretization 
process, one can find a discrete system which has no PN barrier. One would expect kink dynam- 
ics in such a system to be much simpler and more continuum-like than in the conventional discrete 
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system, even on relatively coarse lattices. This expectation has been confirmed for a similar discrete 
sine-Gordon system Q. The rest of the paper is organized as follows: in section |^ we define the 
discrete 0"* system, in section |^ we prove that it has no PN barrier, in section ^ we present numer- 
ical simulations of the system, and compare its performance with that of the conventional discrete 
system, and finally section presents some concluding remarks. 



2 A topological discrete 0^ system 



The continuum (j)'^ model consists of a real valued scalar field : R^^^ M whose time evolution 
is governed by the Lagrangian L = Ek — Ep, the kinetic and potential energy functionals being, 
respectively. 
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To ensure that (j) has finite energy, one imposes that >±las|a;|— >oo. The trivial solutions 
4>{x, t) = ±1 arc the vacua of the model, while a solution interpolating between —1 and 1 (1 and —1) 
is a kink (antikink). Static kinks arc minimals of -Bp[(/)] and can be found by means of a Bogomol'nyi 
argument Namely, if limj,^±oo = il, then 
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So Ep > I in the kink sector, with equality if and only if 
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Equation ^ is called the Bogomol'nyi equation. Note that it is a first order o.d.e., unlike the static 
field equation which is second order. It is easily solved, yielding 



{x) — tanh -{x — b) 



(5) 



where the constant of integration b is identified as the position of the kink centre. 

Upon discretization, x takes values in hZ = {0,±h,±2h, . . .}, where h is the lattice spacing. 
It is convenient to introduce the notation f+{x) = f{x + h) and /-(x) = f{x — h) for forward 
and backward shifted versions of any function /, and A/ = h^^{f-\- — /) for the forward difference 
operator. Ward's idea is to seek a discrete version of the potential energy of the form 
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where D —> d(j)/dx and — > (1 — (fp') in the continuum limit ^ 0, and D and F have the product 



DF = -A — 
3 
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Such a lattice potential energy clearly has the correct continuum limit. The crucial point is that, 
given (|^), the discrete system has a Bogomol'nyi argument exactly analogous to its continuum 
counterpart. Let (f) : hZ — > M have kink boundary conditions. Then 
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using the boundary conditions cf) — s- ±1 as x ±oo. 
only if (j> satisfies the lattice Bogomol'nyi equation, 

The most natural choice of D and F seems to be [Q 



It follows that Ep > ^ with equality if and 
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D = A(j) 



F = + </>+</> + 02), 
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which clearly have the correct continuum limits, and are easily shown to satisfy equation (Q). Note 
that the lattice Bogomol'nyi equation (||) is a first order difference equation, while the static field 
equation, 

is a second order difference equation. We have shown that if a solution of (^) exists with kink 
boundary conditions then it is a global minimal of Ep (within the kink sector) with energy |, and 
hence a static kink solution. The next section is devoted to proving the existence of such solutions. 

We will refer to this lattice system as the topological discrete system (TDt/f^S) since it preserves 
the topological lower bound on kink energy. The conventional discrete (p'^ system also has a potential 
energy of the form (|^), but while D is the same (a standard forward difference), F is chosen to be 
the same as in the continuum model, F — 1 — (fi^. In this case, the Bogomo'nyi argument is lost. So 
the essential difference here is that we have "distributed" the double well substrate potential over 
pairs of nearest neighbour lattice sites. In both discrete systems, one completes the dynamics by 
defining the obvious kinetic energy functional. 
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Other interesting choices can be made B, but this is the simplest. 



3 The kink moduli space 

We start with a definition: a two-sided sequence <j) : hZ ^ M is a static kink if it is a solution of 
the Bogomol'nyi equation (^), is monotonic and converges to 1 at oo and —1 at — oo. The space 
of all static kinks is called the kink moduli space, denoted, for a particular value of h, by 
Given that discretization has broken the continuous translation symmetry of the continuum model 
to symmetry under integer translations by h, one would expect to be discrete. This is indeed 
true for conventional discrete systems, where kinks must lie halfway between or exactly on top of 
lattice sites, the difference in energy between these two types of static solution being the origin of 
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the PN barrier, which so complicates the kink dynamics. However, the TD0^S , remarkably, has a 
continuous moduli space: static kinks may take any position relative to the lattice. All such kinks 
saturate the Bogomol'nyi bound on energy, so there is no PN barrier resisting their propagation 
through the lattice. Specifically, we will prove. 

Proposition 1 Ifh£ (0, 2] then for all (j)Q G (—1, 1) there exists a unique static kink with 0(0) = (j)o. 

Note that a continuous translation orbit of kinks can equally well be parametrized by kink position 
6 G R (defined, for example by linear interpolation) or by the value of the field at cc = 0, 0(0) G 
(—1,1). So proposition 1 could be rephrased: for all h G (0,2], Mh = R. This is the important 
result for our purposes. It is also interesting to consider what happens to the kink moduli space for 
h > 2. In fact, except for some special values of /i G (2, 3], Mh = 0. However, for an infinite discrete 
subset of (2,3], Mh = hX. Precisely, 

Proposition 2 For h > 3 there exist no static kinks. There exists a strictly decreasing sequence hn 
with ho = 3 and lim„^oo hn — 2 such that if h — /i„ there exists exactly one static kink ( modulo 
lattice translations). This kink is symmetric about its central site(s) and has exactly n -\- 1 sites 
different from ±1. 

Proof of proposition 1: Since the Bogomol'nyi equation, 

^ + 0) 0+ + - ^0 - 3) = (13) 

is forward backward symmetric, that is, invariant under (0_|_,0) 1— > (—0,-0^), it suffices to prove 
that for any h G (0, 2] the right hand sequence converges to 1 for all ^(O) G (—1, 1). Note that given 
G (—1, 1), equation ( p^ ) does not uniquely determine 0+ since both roots are real. However, the 
lower root is always less than 0, and we seek monotonic solutions, so we can discard the lower root 
at each iteration. It follows that if a static kink with 0(0) = 0o G (—1; 1) exists, it is unique. 
It remains to prove convergence to 1 of the sequence 

0(n/i) = r(0o) (14) 

generated by iterating upper root of equation ([l3|), starting at any value 0o £ (^Ij 1) ^^'^ assuming 
that < ft, < 2. Here /" denotes the n-th composition of / (so /" ~ Id), and 
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A straightforward calculation shows that on [—1,1], f{q) > q, the only fixed points being ±1. 
Furthermore, f{q) > 1 if and only if g G Ih, this being the interval (7 — I7 — 1|,7 + I7 — 1|) where 
7 = 3h-^ - 2-1. Since h n [-1, 1] = if ft < 2, we find that / is bounded above by 1 on [-1, 1]. So 
the sequence /"(0o) is monotonic and bounded, and hence, by the monotone convergence theorem, 
convergent. Call its hmit L. We have established that / : [—1, 1] [—1, 1], so the sequence /"(0o) 
lies in the interval [—1, 1]. Either there exists N ^ N such that /"(0o) — L for all n > iV or L is an 
accumulation point of [—1,1]. In either case L G [—1,1] since [—1, 1] is closed. Consider then the 
image of the sequence /"(0o) under /. Since / is continuous, /(/"(0o)) — > f{L). But /(/"(0o)) is 
a subsequence of f^{(j>o), so converges to L. Hence L is a fixed point of /, and since L cannot be 
— 1, it follows that (j){nh) ^ 1. □ 



It is interesting to compare the situation with that of the topological discrete sine-Gordon 
(TDSG) system B. There, one is fortunate enough to have an explicit formula for the static kink 
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solution, obtained by slightly modifying the continuum kink. One finds that Mh = M for all /i < 2, 
but the "proof" of this fact rests on the observation that the explicit formula has a continuous real- 
valued parameter in it. Unfortunately, no continuum-inspired ansatz seems to work for the TD(/)*S 
Bogomol'nyi equation, so we are forced to give a less constructive proof. Note that in both systems 
the kink moduli space degenerates at /i = 2. This may seem an artificial coincidence which would 
be lost if different normalizations were used, but in fact the TDi/i'^S was normalized specifically so 
that the minimum phonon frequency (the angular frequency of small amplitude travelling waves) is 
unity, as in the TDSG system. One then finds that both systems have the same phonon dispersion 
relation, 

-=^-^-oskh (16) 

obtained by substituting a travelling wave ansatz cos(fea; — cut) with frequency lo and wavenumber 
k into the linearized equation of motion (linearized about the vacuum (j) — 1, not (/) — 0). So in 
both systems Mh degenerates precisely where the phonon dispersion relation collapses to a flat line. 
Since there is no obvious relation between the existence of kinks, a nonlinear phenomenon, and the 
radiation of the system, a linear phenomenon, this coincidence is rather strange. One difference 
between the two systems is that, whereas there are no kinks for /i > 2 in the TDSG system, there 
are kinks in the TD^^'S for certain discrete values of h in (2, 3], as stated in proposition 2, to whose 
proof we now turn. 



Proof of proposition 2: We adopt the definitions and notation used in the proof of proposition 1. 
As stated above, f{q) > 1 if q G /h, so a static kink must have 0(0) ^ Ih. Consider the backwards 
iteration (/>_ = /(</>). By forward-backward symmetry, f{q) < —1 if q G /;t = (—7 — I7 — 1|, —7 -I- 
I7 — 1|), where 7 was defined above. So any static kink must also have 0(0) ^ Ih- If /i > 3 then 
(— 1, 1) C //i U //j. Hence there exists no static kink for /i > 3. 

Let h G (2,3] and (j){x) be a static kink. Since lima;_,±oo = ±1, either (j>{±x) — ±1 for all 
X > Nh for some N G Z"*", or either 1 or —1 is an accumulation point of the sequence. The latter 
cannot be the case since (p is bounded away from 1 and —1 by the intrusion of Ih and Ih into [—1,1] 
respectively (see figure 1). So the static kink must be constant away from its centre. This is possible 
if and only if there exists n G N such that 

r(-7+l7-l|)=7-|7-l| (17) 

where /" again denotes the n-th composition of /. For then if 0(0) = —7 -I- I7 — 1|, the backward 
sequence jumps immediately to —1, which is a fixed point, and the forward sequence leads, after n 
iterations to 7 — I7 — 1|, which is subsequently mapped to the fixed point 1. Clearly such a kink is 
unique modulo lattice translations, is symmetric and has n+l sites different from ±1. For example, 
if h — 3 then —7 -I- I7 — 1| = so there is a static kink with 0(0) = and (f){x) — x/\x\ elsewhere. 
The only other static kinks are obtained by lattice translation of this one. 

So for each n G N there is a n+l site kink in the system with spacing h — hn where /i„ is a solution 
of equation (p7|); /iq = 3 is a solution with n = 0. We seek to prove the existence of a decreasing 
sequence of such solutions, (ft-n), converging to 2. Let = 7 — I7 — 1|, so /i G [2, 3] =J> x £ [0, 1], and 
define g{x) = — /(— x)j a graph of which is presented in figure 2. Then the condition ( p!7[ ) for the 
existence of a static kink can be rewritten 

5"(X) = -X- (18) 

Note that g(0) = —1, g{l) = 1 and g is continuous, so g has a root, call it ri, between and 1, 
and so there exists a solution xi of g{x) — —X i^i (O,?"!)- Now, if g" has a root r„ G (0,1), then 
(7"+^ has a root r„+i G (r„, 1) since g^^^{rn) — —1 and g"^^(l) — 1. Hence, by induction, for each 
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n G Z^, g" has a root r„ G (?'n-i, 1), and since ^"(rn-i) = — 1 there exists a solution Xn of equation 
( p^ with e (''n-iifn)- The sequence (xn) is monotonic and bounded, hence convergent, and 
hm„^oo Xn = lim„-+oo r„. 

It remains to show that hm„_+oo = 1 , for then the corresponding sequence of h values converges 
to 2 as claimed. Note that g is convex and increasing on [0, 1], whence it follows that 5 < ri < 1 
(consider the cord from (0,-1) to (1,1) in figure 2). Consider now the function 5" : [r„_i,l] 
[— 1, 1] where n > 2. Defining k :— g""^, one finds that 

(5")"(x) = g"{Hx}Kk'{x}r + g'iKxWix)- (19) 

The first term is nonnegative since g is convex on [0,1], and if k is convex on [r„_2,l] then the 
second term is nonnegative since g' o k > on [r„_i,l] (g is increasing on [0,1]). Hence, by 
induction, : [r„_i, 1] [0, 1] is convex, for each n G Z^. It follows that 1 — 2^" < rn < I and 
hence that lim„_,oo I'n — ^- ^ 



For h G (2,3] then, Mh is discrete. One should not, however, interpret this as the appearance 
of a PN potential in the system, since the discreteness stems from our requiring that static kinks 
be monotonic (so that no site may lie outside [—1,1]). Numerical solution of the lattice Bogoml'nyi 
equation suggests that even if h > 2 there are solutions with the correct boundary behaviour (and 
hence Ep = |) for any (pQ G (—1,1), but that they "overshoot" and approach the vacua from the 
wrong sides (i.e. as x — > 00, (j){x) — > 1 from above, rather than below). We do not regard these as 
being real kink solutions since they are not even qualitatively close to continuum kinks sampled on 
the lattice. The conclusion is, then, that to simulate the continuum theory one should use a lattice 
with h < 2. Note that choosing, for example, h — 1.5, the system is in an extremely discrete regime 
since the kink structure is spread over very few lattice sites (see figure 3). 



4 Numerical simulations 



Having established that kinks of the TD0^S experience no PN barrier, the question arises whether the 
resulting kink dynamics is consequently simpler than that customarily observed in discrete systems. 
We address this question by performing numerical simulations of the lattice equation of motion. 
This involves approximately solving the following system of coupled nonlinear o.d.e's. 
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on a large (but finite) lattice, with a Lorentz boosted continuum kink solution as initial data. 
(Although it would be straightforward to generate a genuine discrete static kink for the initial 
profile, there is no obvious way to boost it, without an explicit formula. In any case, the point of 
the endeavour is to simulate continuum kink dynamics.) Simulations were performed with a variety 
of lattice spacings and initial kink velocities. The solutions were generated using a simple fourth 
order Runge-Kutta scheme. 

In every case it was found that the kink moves freely through the lattice, without being pinned 
(see figure 4), exciting small amplitude travelling waves which propagate backwards in its wake. 
To reduce reflexion of this radiation from the fixed left hand boundary during long simulations, 
the first few lattice sites were damped. Figure 5 presents plots of velocity against time for several 
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long duration simulations. (The thickness of the curve is due to velocity oscillations as the kink 
passes from cell to cell. This dynamical wobble is not caused by any PN type potential, but can be 
understood in terms of geodesic motion on the moduli space Mh with respect to a periodic metric 
[Q.) After an initial drop in velocity as the kink relaxes to a profile more suited to the discrete 
system, subsequent kink deceleration is very modest, even at high speed {v = 0.6) and on coarse 
lattices. 

This behaviour should be contrasted with that of the standard discrete system, whose equation 
of motion, 

0=^(0+-20 + 0_) + i0(l-(/.2) (21) 

was solved using the same algorithm and initial data, for purposes of comparison (also figure 5). 
These results are similar to those obtained by Combs and Yip |^ in their work on the conventional 
discrete (p'^ system. One sees that at every value of h, kink dynamics in the TD(/)"'S is far less 
dissipative, the difference being more pronounced for larger h. Indeed, a kink set off with a speed 
V — 0.6 on a lattice with h — 1.8 quickly becomes trapped in the conventional system (figure 6), but 
remains free (apparently) indefinitely in the TD(/)^S . 



5 Concluding remarks 

In this paper we have shown that by choosing a discretization which preserves the Bogomol'nyi 
bound on kink energy, one can find a discrete 0^ system which has no Peierls-Nabarro potential. 
The resulting kink dynamics is much simpler than one expects of a discrete system, even on coarse 
lattices. The kinks are never pinned, and suffer only modest radiative deceleration. So the TD0^S 
is an efficient and natural choice of discretization for numerical study of the continuum (f)"^ model. 

In higher dimensional field theories with topological solitons (the 0(3) sigma, Skyrme and abelian 
Higgs models for example) the reduction in computational cost would be quite significant if simi- 
larly natural discretizations could be found. Here, however, there is an extra concern: solitons on 
high dimensional lattices tend to be unstable because they can unwind and "fall through" a single 
plaquette. One way to preserve the stability of the continuum model is again to preserve the Bo- 
gomol'nyi bound on soliton energy, an approach pursued by Ward 1^ for two planar field theories. 
In this case it seems impossible to eliminate the PN barrier entirely, although one can attempt to 
reduce it as much as possible. It would be interesting to see to what extent these discrete systems 
have continuum-like soliton dynamics, as has been found here for 0* kinks, and elsewhere ^ for 
sine-Gordon kinks. 
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Figure captions 



Figure 1: The shaded region, excluding boundaries, represents values which get mapped outside 
the interval [—1, 1] under either forward or backward iteration of the Bogomol'nyi equation. The 
intersection of any vertical line with this region is I^U 1^. 

Figure 2: The function g{x} of the proof of proposition 2. 

Figure 3: A static kink of the model with h = 1.5. Circles represent the kink profile itself, while 

the solid line represents its energy distribution, which should be regarded as residing in the links 
between lattice sites. Note the extreme discreteness of this structure: over 99% of the kink's energy 
is carried by the 4 central links. 

Figure 4: Propagation of a kink through the lattice of spacing h = 1, with an initial speed v = 0.2. 
Note the apparent lack of radiation in the kink's wake. 

Figure 5: Comparison of the motion of fast kinks (initial velocity v = 0.6) in the TD(/)^S with 
that in the conventional discrete system, for lattices of spacing h = 1.0, h = 1.4 and h = 1.8. In 
each case velocity is plotted against time for both systems on the same graph, the upper curve 
showing the TD(/)'*S data and the lower curve showing the conventional system data. The difference 
is particularly striking for h = 1.8, since the kink quickly becomes pinned in the conventional system, 
while it remains free indefinitely in the TDt^^S . 

Figure 6: Kink pinning in the conventional system with h = 1.8. Set off at a speed of 0.6, the 
kink manages to travel through only 8 lattice cells before being trapped by the PN barrier. This 
should be compared with the dynamics of a TD(/)'*S kink on the h = 1.8 lattice, which remains free 
indefinitely and after 3000 time units has decelerated by less than 17%. 
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